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A problem of great concern in aviation and submarine propulsion 
is the control of the boundary layer and, in particular, the methods 
to extend the laminar region as a means to decrease noise and fuel 
consumption. In this paper we study the flow of air along an air- 
foil when a layer of ionized gas and a longitudinal electric field are 
created in the boundary layer region. By deriving scaling solutions 
and more accurate numerical solutions we discuss the possibility of 
achieving significant boundary layer control for realistic physical pa- 
rameters. Practical design formulas and criteria are obtained. We also 
discuss the perspectives for active control of the laminar-to-turbulent 
transition fluctuations by electromagnetic field modulation. 
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Whether a flow is laminar or turbulent, the effects of the viscosity of the 
fluid are greatest in regions close to solid boundaries. The region close to 



the boundary, in which the velocity varies from zero, relative to the surface, 
up to its full value, is called the boundary layer. The concept of boundary 
layer introduced by Prandtl in 1904, was a most significant advance in fluid 
dynamics, in the sense that it simplified the study by separating the flow in 
two parts: (1) the region where velocity gradients are large enough to produce 
appreciable viscous forces - the boundary layer itself - and (2) the external 
region where viscous forces are negligible compared to other forces. From 
the computational point of view the concept of boundary layer also plays 
a significant role because, rather than having to deal with time-consuming 
general purpose finite-element codes, results of comparable precision may be 
obtained by fast and relatively simple finite-difference implicit algorithms. 

When a fluid flows past a solid body, an airfoil for example, a laminar 
boundary layer develops, in general only for a very small distance near the 
leading edge, followed by a transition to a turbulent boundary layer. Never- 
theless, because near the solid wall velocity fluctuations must die out, below 
the turbulent region there always is a laminar sub-layer which in general is 
very small (of the order of a micrometer). The transition from the laminar 
to the turbulent region is controlled by the local Reynolds number, defined 
in terms of the effective thickness of the boundary layer. It also depends on 
the smoothness of the surface and on the external perturbations. The skin 
friction drag is proportional to the gradient of the longitudinal velocity at 
the solid boundary. Because of the mixing properties of the turbulent layer, 
the gradient in the laminar sub-layer is much greater than the gradient at a 
fully laminar layer. Therefore transition from a laminar to a turbulent layer 
greatly increases the skin friction drag. Another effect to be taken into ac- 
count is the separation of the boundary layer, which occurs at points where 
the pressure gradient along the surface reverses sign. The eddies, generated 
by the resulting reverse fiow, disturb the fiow and form a wake where energy 
dissipation decreases the pressure, thereby increasing the pressure drag. 

Because of the very large ratio between laminar and turbulent skin fric- 
tion drag, much effort has been devoted to develop techniques to delay the 
transition as a means of decreasing fuel consumption and noise. Care should 
however be taken because, in general, a turbulent boundary layer is more 
stable towards separation than a laminar boundary layer. Some of the active 
control techniques that have been proposed include suction of slow-moving 
ffuid through slots or a porous surface, use of compliant walls and wall cool- 
ing (or wall heating for liquids). Injection of fast-moving fiuid, on the other 
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hand, is effective in avoiding separation but it increases turbulence. Most of 
these ideas are fairly old (see for example [||],0) however, in view of their 
interest for the applications and to obtain a more accurate characterization 
of the physical mechanisms, studies of boundary layer control using these 
aerodynamic methods are still, at present, being vigorously pursued (see for 
example || [| § and papers in 0) 

Another class of techniques for active boundary layer control consists in 
acting on the flow by means of electromagnetic forces. Here different tech- 
niques should be envisaged according to whether the fluid is weakly conduct- 
ing (an electrolyte like seawater or an ionized gas) or a good conductor (like a 
liquid metal). Proposals for boundary layer control by electromagnetic forces 
are also relatively old and trace its origin at least to the papers of Gailitis 
and LielausisQ, Tsinober and ShternP] and Moffat [|1^ in the sixties. Inter- 
est in these techniques has revived in recent years and some more accurate 
calculations and experimental verifications have been carried out, mostly in 
the context of electrolyte fluids [|lT| |jl2 . 



In this paper we will be concerned with the flow of air along an airfoil 
when a layer of ionized gas is created on the boundary layer region. Local 
ionization of the air along the airfoil is not practical from the technological 
point of view, therefore we will assume that a stream of ionized air (or some 
other ionized gas) is injected through a backwards facing slot placed slightly 
behind the stagnation point (Fig.l). The body force that we consider to 
be acting in the ionized fluid is a longitudinal (along the flow) electric field 
created by a series of plate electrodes transversal to the flow and placed inside 
the airfoil with the edges on the airfoil surface. 

The emphasis of our study is on finding physically reasonable ranges of 
parameters and analytic approximations that might lead to simple designing 
procedures. For this purpose, before the numerical calculation of Section 
3, we dedicate some time to the study of scaling solutions and analytical 
approximations. 

The provisional conclusions of our study are that it is possible to use 
this technique to control the profile of the boundary layer laminar region. 
With the rates of ionization that are needed and the injection method, it 
is probably unrealistic to expect that the laminar region may be extended 
over all the airfoil in normal (aviation) working conditions. Therefore this 
method should be used in conjunction with methods for control of turbulent 
boundary layers (riblets, large-eddy breakups, additives, etc.) in the rear part 
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of the airfoil. Also the injection of the stream of ionized gas in the leading 
edge may, by increasing the velocity component normal to the airfoil, create 
turbulence. Therefore it seems advisable to have a compensating suction 
region after the injection slot. The ionization rate will also improve if the 
gas extracted through the suction region is recycled through the ionizer. 

Notice also that, once the fluid in the boundary layer is ionized, large 
scale velocity fluctuations may be detected by a few local probes. This raises 
the possibility of obtaining a negative feedback effect by an appropriate time- 
dependent modulation of the electric field. By controlling the growth of the 
velocity fluctuation in the transition region, a further extension of the laminar 
region may be obtained. This is briefly discussed in the last section of the 
paper. 

The overall conclusion is that, when used in conjunction with other tech- 
niques, as explained above, the method of boundary layer control by electric 
flelds might be interesting from the fuel consumption point of view. This 
study was carried out as a preparation for an experiment being set up in our 
Mechatronics Laboratory. 



2 Ionized boundary layers with electric fields 
2.1 The boundary layer equations 

We use orthogonal curvilinear coordinates with x parallel to the surface along 
the flow and y normal to the surface. If k5 is small {k denoting the curva- 
ture and 5 the boundary layer thickness) the conservation and momentum 
equations in the incompressible fluid approximation may be written 

dv ^dv ^dv 1 dp _ 
dt dx dy Pn 
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u and V are the components of the fluid velocity fleld along the x and y 
directions. Pm is the mass density, the electric charge density and E an 
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applied electric field. The tilde denotes quantities in physical dimensions to 
be distinguished from the adimensional quantities defined below. We consider 
typical values L^, 5r, Ur, Pr,i^r,(^r, Er as reference values for, respectively, the 
airfoil width, the boundary layer thickness, the fluid velocity, the fluid mass 
density, the kinematic viscosity, the fluid charge density and the electric field. 
Then we define the adimensional quantities 

X y u vLr , ,s 



Pm P Tt Ur^r /r\ 

Pr PrUr 



u^-,a^-, E=- (6) 



1 S"^ 

In general Rl » 1- Neglecting terms of order and we obtain, for 
stationary solutions = f| = 0^ 

dx dy 

du du 1 dp d^u 1 , , ^ , . 
dx dy pmox dy^ pm 

^ = ^-fcr{x, y)Ey{x, y) (9) 
dy Lj. 
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where uj = = and 7 = h^f^. Unless the electric field component 

normal to the airfoil is very large, one has |^ and the pressure term in 
the second equation may be expressed in terms of the fluid velocity far 
away from the airfoil 

du du due d'^u 1 , ,^ , , 

+ = + ^^^r^ + ^ — yj^A^^ y) (10) 

ox oy ox oy^ pm 

To take into account turbulence effects one should also replace in (|T^) the 
velocity fields u and f by m + m and v + v , u and v being fiuctuation fields 
with zero mean, u' = 0, v' = 0. The effect of the turbulent field on the mean 
flow is now obtained by taking mean values. In a two-dimensional turbulent 
boundary layer the dominant eddy stress is —u'v. Assuming the eddy shear 

du 
dy 

-^-=% (11) 



stress —u'v' and the mean rate of strain |^ to be linearly related 



one obtains flnally 



du du due d ( ^du\ 1 , , ^ , , 

+ = + v^-^ Ptt- + 7 — cr{x, y)E^[x, y) (12) 

ox oy ox oy \ oy J p„ 



-'m 



with 

/? = 1 + - (13) 

being, in general, a function of y through the dependence of the eddy viscosity 
e on the local velocity fleld. (3 should be obtained from a turbulence model. 

To analyze the scaling solutions and for the numerical calculations in Sect. 
3 we deflne a stream function ip and make the following change of variables 

\uujj 4(x) 

^ = {ueUuj)^i{x)f{x,Ti) (15) 

u = ^ v = -— (16) 

dy ' dx 
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The continuity equation is automatically satisfied by (0) and one is left 
with 



drj \ drf I dx drf Ue dx \ drf dx dr] drjdx I u^pm ^ x , 



(17) 



2.2 Scaling solutions 



We assume that the electric field to be created by a series of plate electrodes 
along the z-direction, that is transversal to the fluid flow. For this electrode 
geometry the mean electric field in the x-direction may be parametrized by 

where x and y are the adimensional coordinates defined in (Q). Eq = |M 
is the field at ?/ = 0, controlled by the potential differences between the 
electrodes, and l{x) is of the order of the electrode spacing. For thin (laminar) 
boundary layers the field may with good approximation be considered to 
be independent of y throughout the boundary layer thickness, as long as the 
appropriate charge density profile is chosen (see below). 

For the main application we are addressing, ionized air would be injected 
through a slot near the leading edge of the airfoil, being then carried along 
the airfoil surface by the flow. The steady-state charge distribution in the 
boundary layer is obtained from the continuity equation 

dci da . 

—u + —v = j 19 
ox oy 

j being the source of electric charge. For a point source at the position 
{xo, yo), that is j = c(xo, yo)S{x — xo)6{y — yo), the solution is 

a{x, y) = c(xo, yo)5 (^/'(x, y) - ^'(xo, yo)) 0{x - xq) (20) 

ip being the stream function defined before. Then, for a column of ionized air 
injected at a backwards facing angle through a slot placed at xq, behind the 
stagnation point, each point acts as a point source of intensity proportional 
to the local fluid velocity. Furthermore the intensity of the effective source 
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is depleted up the column. Taking the depletion effect into account, one 
obtains by integration of Eq. (|20|) 

a{x, y) = do (1 - diij{x, y)) 6* (1 - y)) 9{x - xq) (21) 

(Jo is the injection intensity and di characterizes the rate of depletion. The 
conclusion is that the charge density is maximum at the airfoil surface, de- 
creasing to zero at a distance that depends on the fluid dynamics and the 
injection regime. In numerical simulations one may easily use the fairly 
accurate equation ( pT]) for the charge density profile. Here however, the 
dynamically-dependent charge density profile will be parametrized by the 
simpler formula 

a{x,y) = cTo(l--) (22) 

We now look for scaling solutions of (p^. A scaling solution is one for 
which / is only a function of rj. Eq.(^) becomes 

{Pf")\i I sf+^-^ = -^e('^)o. (1 - /) „, , f „, , , 

^ ^ Ue ox UePm ^ ' uj'^{x) + UJiyt,^[X)ri^ 

(23) 

with boundary conditions 

/(0) = /'(0)=0 /'(oo) = l (24) 



where, for simplicity, we have denoted /' = |^ and ^= 

Let the pressure be approximately constant for length scales L of the 
order of the airfoil, that is ^ ~ 0. Let also /? be a constant. This is the 
case for the laminar part of the boundary layer. Then the factorized nature 
of Eq.(^) implies that solutions exist only if 

Ci, Cs and C4 being constants. Therefore 

^(x) = y/CiX + C2 

9i^) = #) (26) 
l(x) = v/cl^(x) 
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There are two physically interesting situations. The one with ci 7^ C2 = 
and the one with ci = C2 7^ 0. The first one corresponds to a boundary 
layer starting at x = and growing with and the second to a constant 
thickness boundary layer. The first one corresponds to an equation 

f"\v) + \f{ri)f"{v) + (1 - f\ri)) = (27) 



with ci = /3 , C2 = , = »V£4 ^ ^j^g second to 

/"'(^) + (l-/'(^))^ = (28) 

with c, = 0,C2^0,a = Sg;g,6=y^. 

In the first case one chooses C2 = to obtain a boundary layer starting at 
a; = 0. The scaling hypothesis requires then an electric field that is singular 
at X = 0, y = {Ex ~ x~^). In any case this electric field solution is not very 
interesting for our purposes because it leads to a boundary layer growth of 
X2 , as in the free force Blasius solution. Therefore it will be more interesting 
to consider a small field free region in the leading edge of the airfoil and match 
the Blasius solution there with the constant thickness solution of Eq.(^). 

Gailitis and LielausisQ have also obtained a theoretical solution of con- 
stant thickness. However they consider a different force field distribution and 
no dependence of the fluid charge density on the boundary layer dynamics. 
Therefore their boundary layer profile has a very different behavior. 

The solution of Eq.(|28|) is easily obtained by numerical integration (see 
below). Notice however that with the replacement 

0(r/) = l-/'(r/) (29) 



and choosing 04 = ^, which is a simple rescaling of ^, Eq .(|28D becomes the 
zero-eigenvalue problem for a Schrodinger equation in the potential a/{l+T]'^), 

-0"(^) + T^0(^) = O (30) 
1 + r]^ 

One may use the well-known WKB approximation to obtain 

/'(") - 1 - , (31) 

r] + ^/l + Tf 
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Eq.(|3T|) is a very good approximation to the exact solution for a > 1 (see 
Fig. 2). Fig.3 shows the effective boundary layer thickness as a function of 
a. The effective boundary layer thickness 5* is defined here as the value of r] 
at which the velocity u reaches 0.95 of its asymptotic value u^- A very fast 
thinning of the boundary layer is obtained (several orders of magnitude) for 
a relatively short range of the a parameter. Fig.3 shows the variation of 5* 
for small a. For large a (and small 5*) one has the asymptotic formula 

„ 2.9957 

which is obtained from Eg . (|3ll) . 

If the longitudinal electric field is assumed to be a constant {Eq) 
throughout the boundary layer thickness, with the same charge profile, the 
solution is even simpler, namely 

/'(r^) = 1 -e"''^ (32) 

with = ^/c2 and 

Again, since ^ is a constant, this is not fully realistic because it leads to a 
constant thickness boundary layer. 

For reference values of the physical quantities in Eqs.(H-|3) we take 



Ur = 100 m s"^ 
Lj. = 1 m 
6r = 10-3 m 

= 1.2 Kg m-3 (34) 

Er = 500 V cm"^ 
(T.r = 15 fiC cm^3 
Ur = 1.5 X 10"^ m^ textnormals~^ 

For these reference values, the adimensional constants uj and 7 defined after 
Eq.d) are 

a; = 0.15 

7 = 62.499 ^ ^ 
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For comparison we mention that in the classical force-free Blasius solu- 
tion, and for these reference parameters, the y— coordinate y* corresponding 
to 5* (that is, the point at which ^ = 0.95) is 

y* = 1.55 X lO^^v^ (36) 

Stability of a laminar boundary layer cannot safely be guaranteed for local 
Reynold numbers greater than about 10^. Therefore requiring 

= ^ ^ 10' (37) 
u 

one obtains, for the reference parameters, y* ~ 0.15 mm. Using (|3^) the 
conclusion is that, for these parameters, the laminar part of a force-free 
boundary layer is only of the order of 1 cm, just a tiny portion of a typical 
wing. 

Now we use the scaling solutions ( p2D and (|3TD to obtain an estimate of 
the effects of a longitudinal electric field. For the constant field case (^) 
from 

f\S*) = 1-6-^*^ = 0.95 

and 

2/* = <5*v^^/^ = 0.15 

V We 

using (^) one obtains 

do = 0.957 

That is, to insure a constant thickness boundary layer with local Reynolds 
number Rs = 10^ (at the point where ^ = 0.95), one needs a charge density 
ao at y = 0, in physical units (and for the reference values of the kinematical 
parameters) 

ao = (ToCr = 14.36 /iC cm~' 

For the variable field case (pi]) the estimate depends on the separation of 
the electrodes. Taking l{x) = 10, that is an electrode separation of the order 
of one centimeter, and the references values for all quantities except for the 
charge density (namely Eq = = 1, Ue = 1, etc.) one obtains C4 = 0.15, 
g{x) = 10, C3 = g{x)$,{x) = 258.2, = ^{x) = 25.8, and requiring 




= 0.15 



fin = 0.95 
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one finally obtains a = 39887.77 leading to 



(To = 0.957 

the same estimate as above. The large value of a that is obtained shows 
that the WKB expression (|31|) is a good approximation for physically inter- 
esting parameter values. On the other hand the fact that the same charge 
density estimate is obtained both in the constant-field and the variable-field 
cases, shows that it is realistic to consider the field as approximately con- 
stant throughout the laminar boundary layer thickness, as long as a variable 
charge profile (|2l| ) or (|2^) is used. 



The above estimates were obtained using the reference values for the 
kinematic variables. For other values we have the following designing formula 
(in normalized units) 

3 Numerical results 

For the numerical solution of Eq.(|l^), with a given by Eg. (^2]), we use an 
implicit finite-difference technique ([|1^ - ||16[). Define F{x,ri) by 



F{x,v) = ^ (39) 



and 



(40) 



^ uipmP^ " UefJ ax 

Then Eq.(|T^ becomes 

9''F dF ^ dF ^ 

TTT + + ^2^^ + asT^ + ^4 = (41) 

orj'^ or] ox 

The derivatives are replaced by finite-difference quotients with a variable grid 
spacing concentrated near rj = 0, where F changes more rapidly. Let k > 1 
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be the ratio between two successive grid spacings in the 77— direction. 

k 



Vi+i - Vi 



Then 



(dF\ 
\dv Ji+1 



+1J ^ ^2 

Ai 



\dxJi^lj Ax 
Al = r]j+i - Vj + k^iVj - Vj-l) 



Substitution in Eq.(^T]) yields 

+ BjFi+ij + DjFi+ij^i + Gj = (42) 



with 



A . — _2 I ai 

- A2 + Al 
_ -2(l+fc) _ ai(l-fc2) , „ , 03. 
~ A2 Al "2 

2k aik"^ 



D = — 



The boundary conditions aX rj = and 77 — 00 are known 

/(^,l)=F(z,l) = 
F{i,N) = 1 

where is the largest label of the grid, in the r/— coordinate, chosen to be 
sufficiently large. 

Because of the tridiagonal nature of the solution in the line i + 1 is 
obtained by the two-sweep method, the recursion relations being 



with ai = and f3i = 0. 

To start the integration process there are basically two methods. In the 
ffist the integration is performed from left to right in the a;— coordinate with 
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the grid extended to the left of the airfoil, where the flow is known. With 
the solution known in the line i, the coefficients Aj to Gj for Eg . (^21) are 
computed at the point Notice that f{i,j) is obtained by integration of 

the solution F. ^ 

f{t,r^)= [ Fit,C)dC 



The integration now proceeds along the lines, from left to right. After a 
complete pass the process is restarted using now for the calculation of the 
coefficients Aj, Bj, Dj and Gj the old values of F at (z + 1, j). The process 
is repeated several times until the solution stabilizes. 

In the second method, which is the one we actually use, the integration 
process starts from an approximate solution. The scaling solutions derived 
in Sect. 2 are particularly useful for this purpose. 

For our calculations we considered an electric field parametrized as in 
Eq.(IT8|), namely 



-Err = El 



with = 666.66 which corresponds to / = 10, Me = 1 and vuj = 0.15. 
Notice that for these parameters, as pointed out in Sect. 2, the electric field 
has only a small variation throughout the boundary layer region. For the 
scaling function we take ^(x) = ^/x and consider /5 = 1. Then all results 
depend only on the variable S 

1 7 
62.499 ujpm 

{S = 1 when all quantities take the reference values). 

In Fig. 4 we show a contour plot of the numerical solution for / (x, 77) 
(=— ) when S = 0.6. From the x— dependence of the numerical solutions we 
may compute the effect of the electric field in extending the laminar part of 
the boundary layer. By defining, as in Sect. 2, the length of the laminar part 
as the X— coordinate corresponding to a local Reynolds number of 10^ and 
denoting by Xq (;;^ = 0.95) the force-free value we have obtained for the ratio 

Xo 

the results shown in Fig. 5. For = we obtain the Blasius solution and 
as we approach 5* = 0.957, corresponding to the scaling solution, the ratio 
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diverges. The matching of the results in the force-free and scahng hmits is a 
good check of the numerical algorithm. A clear indication of the results in 
Fig. 5 is that not much improvement is obtained unless one is able to obtain 
ionization charge densities of the order of the reference value a^- 

4 Discussion and conclusions 

^ In this paper we have concentrated on controlling the profile of the bound- 
ary layer. The profile has a direct effect on the laminar or turbulent nature 
of the flow which, in a simplified manner, we estimated by a local Reynolds 
number (^) defined as a function of the effective thickness. Another relevant 
aspect, of course, is the active control of the transition instabilities that can 
be achieved by electromagnetic body forces on the charged fluid. 

Turbulence and transition to turbulence are three-dimensional phenom- 
ena. However, for the large scale small amplitude (Tollmien-Schlichting) 
fluctuations, that first appear in the transition region, a two-dimensional 
model is a reasonable approximation. In Eqs.(|I]-^ we make as before the 
change of variables , neglect terms of order , ^ and and split the 
velocity and electric fields into 

u = u + u 

V = V + V 

Ex = E,j. + 

where u , v , are the steady-state solutions and u , v' , E'^ the time- 
dependent components. Because of the continuity equation 

du dv' ^ 
dx dy 

we may define a fiuctuation stream function x 

u =— V = - — 
dy ' dx 



Now we assume the fiuctuation to be a (small-amplitude) wave-like function 
of X, y and t 

x(x,y,t)=F(|/)e^(---^*) 
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The imaginary parts of 6 and a control, respectively, the growth rates of 
temporal and spatial fluctuations. The (modulation) electric field is assumed 
to have a similar form 

One now obtains 

/ _ du \dF du _d^F d^F 7 

lua + lO] ^ ta—F + v^—r = z^c^tt^ 

\ dx J dy dy dy^ dy^ pm 

The conclusion is that a space-time modulation of the electric field, with 
the appropriate phase, is equivalent to an effective viscous damping effect 
which delays the growth of the transition region instability. For this to be 
effective one needs to detect the phase of the wave instabilities by electromag- 
netic probes. Absolute synchronization of the feedback electric modulation 
is however not so critical as in acoustic noise cancelation, because here the 
objective is only to obtain an effective damping effect. The simplified treat- 
ment of the transition instabilities is justified by the fact that it is only for 
the small amplitude large scale fluctuations that one may hope to be able to 
detect the phase with some reasonable accuracy. 

^ The kinematic reference parameters defined in (|3^ correspond to typ- 
ical aviation conditions. The conclusion, both from the scaling solutions in 
Sect. 2 and the numerical results in Sect. 3, is that, to obtain a significant 
controlling effect on the boundary layer by this method, the charge den- 
sity (To (at 2/ = 0) should be of the order of the reference charge density 
((Tr = 15yuC cm~^). This charge density corresponds to about 50 times the 
ion concentration a few centimeters away from the emitter of a commercial 
table-top negative corona discharge air purifier with a power of less than 6 
watts. Therefore, it seems technically feasible to achieve a significant bound- 
ary layer control by this method. Another possibility would be to use, instead 
of air, some other easier to ionize gas. This could then be partially recovered 
and recycled by suction. 

7^ As explained in the introduction and because of the perturbation in- 
duced by the injection method, it seems advisable to use this method in 
conjunction with suction and passive control in the rear part of the airfoil. 
Even if a fully laminar boundary layer may never be completely achieved, 
just remember that any small improvement becomes, in the long run, quite 
significant in terms of fuel consumption. 
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7^ The formula (|38|), derived from the scahng solutions, provides rough 
design estimates. Better control over design parameters we hope to obtain 
from the experimental work. 

5 Figure captions 

Fig.l Airfoil transversal cut showing ionized air injection, suction pump and 
plate electrodes. 

Fig. 2 Exact (— ) and approximate (■ ■ ■) constant thickness scaling solution 

f'iv)- 

Fig. 3 Effective boundary layer thickness 6* = ^ = 0.95^ for the 

constant thickness scahng solution. 

Fig. 4 Contour plot of / {x,ri) for S = 0.6. 

Fig. 5 Ratio of boundary layer laminar regions with and without electric 
field control. 
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